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ABSTRACT 

We have made a comparative study of morphological evolution in simulated DM halos 
and X-ray brightness distribution, and in optical clusters. Samples of simulated clusters 
include star formation with supernovae feedback, radiative cooling, and simulation in 
the adiabatic limit at three different redshifts, z = 0.0,0.10, and 0.25. The optical 
sample contains 208 AGO clusters within redshift, z < 0.25. Cluster morphology, 
within 0.5 and 1.0 h“^ Mpc from cluster center, is quantified by multiplicity and 
ellipticity. 

We find that the distribution of the dark matter halos in the adiabatic simulation 
appear to be more elongated than the galaxy clusters. Radiative cooling brings halo 
shapes in excellent agreement with observed clusters, however, cooling along with 
feedback mechanism make the halos more flattened. 

Our results indicate relatively stronger structural evolution and more clumpy dis¬ 
tributions in observed clusters than in the structure of simulated clusters, and slower 
increase in simulated cluster shapes compared to those in the observed one. 

Within z < 0.1, we notice an interesting agreement in the shapes of clusters ob¬ 
tained from the cooling simulations and observation. We also notice that the different 
samples of observed clusters differ significantly in morphological evolution with red- 
shift. We highlight a few possibilities responsible for the discrepancy in morphological 
evolution of simulated and observed clusters. 

Key words: clusters: morphology - clusters: evolution - clusters: structure - clusters: 
statistics 


1 INTRODUCTION 

The hierarchical clustering is the most popular model for the 
Large Scale Structure (LSS) formation. The model relies on 
the assumption that larger structures result from the merg¬ 
ing of smaller sub-clumps. Theoretical paradigm of the hier¬ 
archical evolution is the Cold Dark Matter (CDM) scenario 
which assumes that baryonic matter (stars, hot X-ray gas) 
evolves in the dark matter (DM) potential through violent 
processes. Structural evolution in cosmological objects such 
as galaxies or clusters of galaxies is the underlying principle 
in this scenario. A generic prediction of the CDM model is 
the nonsphericity of the DM halos. The degree of flattening 
of the halos evolves in cosmological time, from highly irreg¬ 
ular at the distant past towards more regular at the present. 
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In principle, the model prediction can be tested comparing 
the DM halo shapes with that of the (baryonic) matter dis¬ 
tributions. A comparative morphological analysis between 
model and observation could help constraining the nature 
of the DM and its role on the LSS. 

Melott, Chambers & Miller (2001; hereafter MCM) 
have reported evolution in the gross morphology of galaxy- 
clusters (quantified by ellipticity) for a variety of optical 
and X-ray samples for z < 0.1. They infer that the evidence 
is consistent with a low matter density universe. Using a 
similar shape measure as well as intra-cluster medium tem¬ 
perature and X-ray luminosity, Plionis (2002) has presented 
evidence for recent evolution in optical and X-ray cluster of 
galaxies for z < 0.18. In both studies evolution is quantified 
by the change of cluster ellipticity with redshift. In a recent 
study, Jeltema et al. (2005) have reported structural evo¬ 
lution of clusters with redshift where cluster morphology is 
quantified by the power ratio method (Buote & Tsai 1995). 
Jeltema et al. used a sample of 40 X-ray clusters over the 
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redshift range ~ 0.1 — 0.8 obtained from Chandra Obser¬ 
vatory. In spite of methodological differences, the results of 
these studies indicate evolution in the morphology of the 
largest gravitationally bound systems over a wide range of 
look-back time. 

The observational evidence prompted concerns about 
the formation and evolution of structures in the CDM sce¬ 
nario via numerical simulations. If the results of simulations 
provide faithful representations of the evolutionary history 
of cosmological objects, then one would expect a similar 
trend in the structure of simulated objects. So far almost 
all studies of simulated clusters are focused either on under¬ 
standing the nature of the background cosmology within 
which the present universe is evolving (Jing et al. 1995; 
Crone, Evrard & Richstone 1996; Buote & Xu 1997; Thomas 
et al. 1998; Valdarnini, Ghizzardi & Bonometto 1999; Suwa 
et al. 2003) or on understanding the distribution and shape 
of the DM halos in various types of simulations, e.g., simu¬ 
lations with or without baryons and gas physics (Dubinski 
& Carlberg 1991; Dubinski 1994; Aninos & Norman 1996; 
Tissera & Dominguez-Tenreiro 1998; Bullock 2001; Buote et 
al. 2002; Jung & Suto 2002; Gao et al. 2004a,b; Kazantzidis 
et al. 2004; Springel, White & Hernquist 2004; Allgood et al. 
2005; Nagai & Kravtsov 2005; Flores et al. 2005; Libeskind 
et al. 2005; Maccio et al. 2005; van der Bosch et al. 2005; 
Zentner et al. 2005). Until recently a comparative study 
of morphological evolution in simulated and real clusters 
was absent. Floor et al. (2003) and Floor, Melott & Motl 
(2004; hereafter FMM) have investigated evolution in clus¬ 
ter morphology simulated with different initial conditions, 
background cosmology, and different physics (e.g. simulation 
with or without radiative cooling). They have used eccen¬ 
tricity as a probe to quantify evolution. Their studies, em¬ 
phasizing shape in the outer regions of clusters, suggest slow 
evolution in simulated cluster shapes compared to the ob¬ 
served one. However, the studies of Floor and collaborators 
are indirect in a sense that they did not analyze observed 
clusters using the same measurement technique applied to 
their simulated data sets. 

In this paper we make a comparative analysis between 
simulated and observed clusters where both data sets are 
juxtaposed and analyzed using the same set of structural 
measures. We analyze cluster morphology and its evolution 
using shape measures such as multiplicity (M) and elliptic- 
ity (e) derived from the Minkowski functionals (Rahman & 
Shandarin 2003, 2004, hereafter RS03 and RS04; Rahman et 
al. 2004). The MFs provide a non-parametric description of 
the images with no prior assumptions made on the shapes 
of the images. The measurements based on the MFs ap¬ 
pear to be robust and numerically efficient when applied to 
various cosmological studies, e.g., galaxies, galaxy-clusters, 
CMB maps etc. (Mecke, Buchert & Wagner 1994; Schmalz- 
ing et al. 1999; Beisbart 2000; Beisbart, Buchert & Wagner 
2001; Beisbart, Valdarnini & Buchert 2001; Kerscher et al. 
2001a, 2001b; Shandarin, Sheth & Sahni 2004). Various mea¬ 
sures, constructed from the two-dimensional scalar, vector, 
and several tensor MFs have been described and tested in 
RS03 and RS04. To derive the parameters applied in this 
study we use the extended version of the numerical code 
developed in RS03 and RS04. 

We study evolution in the simulated clusters in a flat 
CDM universe (ACDM; = 0.3, Oa = 0.7) obtained from 


three different sets of high resolution simulations (Motl et 
al. 2004). The first set has clusters simulated in the adia¬ 
batic limit, the second set contains clusters with radiative 
cooling (RC), and the last set includes clusters with cool¬ 
ing -I- star formation and supernovae feedback (SFF). Fach 
sample contain DM as well as X-ray brightness distributions 
at three different redshifts, « = 0.0, 0.1, and 0.25,. For com- 
parision we also analyze a sample of AGO (Abell, Corwin 
& Olowin 1989) clusters within 2 < 0.25. The sample con¬ 
tains 208 optical clusters derived from 10-inch photographic 
plates taken with the 48-inch Palomar Schmidt Telescope 
(Trevese et al. 1992; Flin et al. 1995; Trevese et al. 1997; 
Flin et al. 2000). 

The objective of our study is twofold: first, to check the 
efficiency of the parameters differentiating various sets of ob¬ 
jects, and second, to explore (statistical) correspondence in 
the morphological properties of the distributions of DM ha¬ 
los, X-ray emitting gas, and optical clusters using measures 
that are sensitive to shape and sub-structures. 

In the CDM model (satellite) galaxies are associated 
with the DM sub-halos that are accreted by their (current) 
parent halo, a bigger structure usually associated with a 
galaxy cluster. If this is the case statistical properties of 
galaxies regarding mass, sub-structure, shape etc., would 
show a similar trend to that of the sub-halos. On the other 
hand X-ray emitting hot gas, evolving in the DM back¬ 
ground potential, would not directly follow the DM distribu¬ 
tion because of its isotropic pressure support. Therefore, a 
statistical analysis of various properties of DM halos, galaxy 
clusters, and X-ray gas distributions will be useful to probe 
possible bias of luminous galaxies toward sub-halos and their 
correspondence with the distribution of hot gas. This is the 
motivation behind the second objective. 

The organization of the paper is as follows: simulation 
technique and the observational data are described briefly 
in §2, a short discussion of shape measures is given in §3. 
The results are presented in §4 and the conclusions are sum¬ 
marized in §5. 

2 DATA 

2.1 Numerical Simulations 

We have analyzed images of simulated clusters projected 
along three orthogonal axes. The clusters have been simu¬ 
lated in the standard, flat cold DM universe (ACDM) with 
the following parameters: fib = 0.026, = 0.3, Ua = 0.7, 

h = 0.7, and as, = 0.928. For a complete description of the 
simulations see Motl et al. (2004). We have used three sam¬ 
ples of clusters derived from the same initial conditions and 
background cosmology. The difference between the samples 
is in the energy loss mechanism experienced by the bary- 
onic fluids. In the first sample, no energy lose is allowed; in 
the second sample fluid is allowed to lose energy via radia¬ 
tion and subsequently cool; in the third sample, physics of 
star formation and supernovae feedback are incorporated in 
addition to radiative cooling. 

The simulations use a coupled N-body Fulerian hydro¬ 
dynamics code (Norman & Bryan 1999; Bryan, Abel & Nor¬ 
man 2000) where the dark matter particles are evolved by 
an adaptive particle-mesh, N-body code. The PPM scheme 
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(Colella & Woodward 1984) is used to treat the fluid com¬ 
ponent on a comoving grid. An adaptive mesh refinement 
(AMR) is employed to concentrate the numerical resolution 
on the collapsed structures that form naturally in cosmo¬ 
logical simulations. The DM particles exist on the coarsest 
three grids; each sub-grid having twice the spatial resolu¬ 
tion in each dimension and eight times the mass resolution 
relative to its parent grid. At the finest level, each particle 
has a mass of 9 x 10® h“^ M©. A second order accurate 
TSC interpolation is used for the adaptive particle mesh al¬ 
gorithm. Up to seven levels of refinement are utilized for 
the fluid component, yielding a peak resolution of 15.6 h“^ 
kpc within the simulation box with sides of length 256 h“^ 
Mpc at the present epoch. Clusters are selected using the 
HOP algorithm (Eisenstein & Hut, 1998) with an overden¬ 
sity threshold of 160. 

A tabulated cooling curve (Westbury & Henriksen 1992) 
for a plasma of fixed, 0.3 solar abundance has been used to 
determine the energy loss to radiation. Heat transport by 
conduction is neglected in the present simulations since it 
has been shown that even a weak, ordered magnetic field 
can reduce conduction by two to three orders of magnitude 
from the Spitzer value (Chandran & Cowley 1998). However, 
Narayan & Medvedev (2001) has shown that if the chaotic 
magnetic field fluctuations extends over a sufficiently large 
length scales within the intra-cluster medium (ICM), then 
thermal conductivity becomes significant to the global en¬ 
ergy balance of the ICM. Energy input into the fluid from 
AGN is also neglected in the current simulations. 

The prescription of Cen & Ostriker (1992) has been 
used to transform collapsing and rapidly cooling gas into 
collisionless star particles. At the finest resolution level, a 
grid cell is eligible to form a star in a given time step if 
the local flow is converging, the dynamical time exceeds the 
cooling time and a Jean’s mass worth of gas exists within the 
cell.To model a population of prompt supernovae, thermal 
feedback has been introduced. The amount of feedback has 
been set from numerical experiments to provide a reasonable 
amount of mass in star particles. The feedback is approxi¬ 
mately 7 X 10^® ergs per solar mass of stars formed or about 
half a keV of energy per particle in the final clusters. 

We have 41 three dimensional clusters from each sam¬ 
ple, giving a total of 123 projected clusters in the respective 
samples. Each projection is constructed within a 8 h“^ Mpc 
(comoving) frame containing 360 x 360 pixels. Majority of 
clusters in each sample is in the mass range ~ 10 ^® — 10 ^“^ 
M© with few clusters (~ 15) in the limit ~ 10^® M©. 

2.2 Optical Clusters 

The details of data acquisition and processing of the opti¬ 
cal sample has been described in Trevese et al. 1992; Elin 
et al. 1995; Trevese et al. 1997; Flin et al. 2000. Here we 
highlight only the essential features of the sample needed 
for this study. 

The sample contains 208 optical clusters, within a < 
0.25, derived from 10-inch photographic plates taken with 
the 48-inch Palomar Schmidt Telescope. It contains rich and 
massive AGO clusters with richness R > 1 and mass, ap¬ 
proximately, in the range ~ 10^® — 10^^ M©. Highly massive 
structures, e.g. Goma cluster (A1654) or clusters constitut¬ 
ing Shapley condensation, are absent in this sample. 


The visual control is the greatest advantage of this sam¬ 
ple. The essential difference between this and other samples 
is visual control of all objects classified as galaxies when au¬ 
tomatic procedure were applied. The visual inspection was 
made for objects with magnitude range at least m 3 -|- 3 mag. 
The relationship between the number of objects with re¬ 
spect to the magnitude and the luminosity function for each 
separate cluster show that clusters are complete at least in 
the magnitude range m 3 to m 2 . 5 . In majority of cases, it is 
complete till m 3 -I- 3. 


3 MORPHOLOGICAL PARAMETERS 


We use multiplicity (M) and ellipticity (e) as quantitative 
measures to study evolution of observed and simulated clus¬ 
ters. Ellipticity is derived from the area tensor functional, a 
member of the hierarchical set of the MFs. This functional 
is given by. 


Aij — 



- M){xj 


Aj) da, 


K 


( 1 ) 


where K is the region bounded by a given contour and Ai 
is the area vector functional, i.e. area centroid, expressed as 
follows, 


Ai = 


As 


Xi da. 


( 2 ) 


The symbol A© represents the area within the contour. It is 
known as the scalar area functional and is given by, 



K 


The area vector functional is in fact the center of mass of 
the region within the contour if we assume that the surface 
density of the (enclosed) region is constant. The tensor Aij is 
closely related to the inertia tensor of a homogeneous region, 
The details of the MFs can be found in Schmalzing (1999), 
Beisbart (2000), and RS03. 

• Multiplicity (M): This parameter is defined as. 


M = 


Aj] 


Y^a. 


As 


(4) 


where Ai is area of the individual components at a given 
level, Amax is the area of the largest component at that 
level , N is the total number of components, and As is the 
total area at that level obtained after summing the areas 
of the components. Multiplicity, M = M{As), is a measure 
with fractional value and gives the number of components 
measured at any brightness level: M = 1 for a single iso¬ 
intensity contour, i.e. component, and M > 1 for multi¬ 
contours. 

It may be mentioned here that Thomas et al. (1998) 
have used multiplicity as a parameter for sub-structure mea¬ 
sure in N-body simulations. They define it as a ratio of mass 
of sub-clumps to cluster mass. In this study it is a ratio of 
the areas (sizes) as defined in equation 1 . 

We use two variants of M to present our results: one is 
the average of multiplicity over all density/brightness levels, 
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Meff, and the other is the maximum of the multiplicity 
found at one of the levels, Mmax- 

• Ellipticity (e): We adopt the definition of ellipticity, 

e = 1 — b/a, (5) 


where a and b are the semi-axes of an ellipse. For our pur¬ 
pose the semi-axes correspond to the “auxiliary ellipse” con¬ 
structed from the eigenvalues of the area tensor (see RS03 
for detail). Note that the “auxiliary ellipse” is an ellipse 
having exactly the same area tensor. 

We have used two variants of e: one is sensitive to the 
shape of the individual cluster components present at a 
given level while the other is sensitive to the collective shape 
formed by all the components present at that level. We label 
these two variants of e, respectively, as the effective (ce//) 
and the aggregate (tagg) ellipticity. Morphological proper¬ 
ties of clusters such as shape and the nature or the degree 
of irregularity existing in these systems can be probed effec¬ 
tively with these two parameters. 

At any given density/brightness level, we construct Ce// 
as a weighted mean normalized by the multiplicity and area 
of the largest contour. 


fe// 


M ■ 


1 

A 


max 


N 

Ai, 

i = l 


( 6 ) 


where ti are ellipticities of the individual components mea¬ 
sured as stated earlier and M is the multiplicity at that level. 
The symbols Ai and Amax have similar meanings as before. 
This measure can be used as an effective tool to quantify 
shapes of large scale merger remnants. 

To construct tagg, we take the union of all components 
present at a given level and form a collective region. The 
integrated region can be expressed as 


R = i?i U i?2 U • • ■ U Rn, 


( 7 ) 


where Ri is the region enclosed by each contour. Subse¬ 
quently we find the components of the area tensor and the 
“auxiliary ellipse” for the region R. 

The behavior of eagg is similar to the conventional ellip¬ 
ticity measure based on the inertia tensor (Carter & Metcalf 
1980). But the construction procedure of these two measures 
are different. The conventional method finds the eigenvalue 
of the inertia tensor for an annular region enclosing mass 
density or surface brightness. On the other hand, the method 
based on MFs finds the eigenvalues of regions enclosed by 
the contour(s) where the regions are assumed to be homo¬ 
geneous. 

We have computed ellipticities after averaging the es¬ 
timates at all density/brightness levels. Our final result is, 
therefore, expressed as Ce// and tagg- 


3.1 Toy Models 

To get a better feeling of the parameters mentioned above we 
provide an illustrative example with toy models. We find this 
demonstration useful since it gives a visual expression how 
the number of group members forming complex structures 
affects the shapes (see also Paz et al. 2005). 

One can think of these toy images as snap shots of differ¬ 
ent clusters (in projection) taken at one particular time. We 
include clusters with different types of internal structures in 


Fig. 1: unimodal elliptic structure (panel 1), asymmetric and 
symmetric bimodal clusters (panels 2 and 3, respectively), 
cluster with filamentary structure (panel 4) etc. The multi¬ 
modal clusters have clumps with different peak brightness. 
We show contour plots of toy clusters at different brightness 
levels where the levels are chosen arbitrarily. For all clusters 
the outer line represents the percolation level where the sub¬ 
structures merge with one another and form a single, large 
system. 

Multiplicity as a function of component area (in grid 
units) is shown in Fig. 2 for our selected toy models. As men¬ 
tioned earlier M is sensitive to the size of the sub-structures. 
The simplest case to see this is a bimodal cluster. For a 
bimodal structure with unequal sub-clumps (panel 2 ), the 
fractional value of multiplicity (1 < M < 2 ) tells us that 
the components of the system have different sizes. The iso¬ 
lated components eventually percolate giving M = 1 at low 
brightness level, i.e, at larger area. On the other hand, for a 
cluster with equal components M — 2 until percolation oc¬ 
curs (panel 3). For clusters with three components (panels 
4 and 5), we see that for a small range of brightness levels, 
the components are well separated where two of these are 
bigger then the third one (2 < M < 3). Afterward two of 
the three clumps merge together giving 1 < M < 2. These 
two remaining components eventually percolate to become a 
single system. The clumps in panel 6 are distributed around 
the center. For this cluster, we see two unequal but well sep¬ 
arated clumps (1 < M < 2) with the same peak brightness. 
The behavior of clusters in panels 7 and 8 is similar except 
that they have a different number of sub-structures. The 
cluster in panel 9 has the largest number of components (a 
total of 7). Two of its clumps are so large compared to the 
other ones that they dominant. The multiplicity is always 
in the range 1 < M < 3 reflecting the merger of clumps at 
different levels. 

Ellipticity for these toy clusters is shown in Fig. 3. In 
this figure the solid and dotted line represent, respectively, 
tagg and teff- For the unimodal cluster in panel 1, teff = 
tagg- For the bimodal cluster in panel 2, the estimate of 
teff is weighted more by the larger component. It is zero for 
the case shown in panel 2. This is also true for the cluster 
in panel 3. However, for a bimodal system with equal sized 
sub-clumps but different elongation, tef / will give an average 
elongation of the two. For systems with sub-structures the 
estimate of tagg, on the other hand, tells us about the overall 
shape of these systems. Due to the presence of two isolated 
components, the system itself appears more elongated than 
the shape of its sub-clumps. An important point to note 
that the estimate provided by tagg depends not only on the 
relative sizes of the components but also on their relative 
separations. This is reflected in the panels containing multi¬ 
clump clusters. For equal separation, a bimodal cluster with 
components similar in shape but unequal in size has lower 
tagg than that of a a bimodal cluster with identical shape 
and size (see the region 1.8 < logj^Q As < 2.2 in panels 2 
and 3 in Fig. 3). In general, as the density and brightness 
level decreases the clumps get bigger and appear closer to 
one another and tagg gets smaller. 

Note that for a multicomponent system with filamen¬ 
tary structure, teff < tagg (panel 4). On the other hand 
if components are distributed around the cluster center, 
Se// > tagg (panel 6 ). The cluster in panel 5 has the unique 
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property that is shown separately by clusters in panels 4 
and 6. In transition at a lower brightness level, the clus¬ 
ter changes its filamentary shape to an extended structure 
where the components are distributed over a region around 
the center. The Cagg profile in panel 8 shows that in the 
range, 2.2 < logj^Q As < 2.8, the cluster develops two, al¬ 
most equal size clumps that are very close to each other. The 
cluster in panel 7 follows the behavior of a bimodal cluster 
except that there is jump in between 2.6 < logj^Q As < 2.8 
where the cluster changes its structure having two unequal 
size clumps to two equal size clumps. The shape of the clus¬ 
ter in panel 9 changes consistently following the merger of 
its clumps at different brightness levels. 

In Fig.s 2 and 3 we also show the variants of the pa¬ 
rameters used later in this study. We use circle and star to 
represent the structural parameters. Me// and Mmax- The 
respective symbols are also used for the shape parameters 
te// and tagg- Comparing this set of parameters with more 
general M and e one can easily see how these measures re¬ 
sponse to the alignment of sub-structures and their spatial 
locations (filamentary, extended etc.). 

3.2 Example of Simulated Clusters 

We demonstrate the behaviors of M and the variants of e 
as a function of area for a collection of simulated clusters 
in Figs. 4 and 5. For each sample we choose two clusters 
at each redshift. We use dark, gray, and faint solid lines to 
represent respectively, the adiabatic, RC, and SFF samples. 
The DM halos and X-ray clusters are shown on the left and 
right panels, respectively. 

Figure 4 shows that both matter and X-ray clusters with 
cooling, generally, have a higher number of sub-clumps than 
those without cooling. Figure 5 shows that in most cases the 
central part of cluster consists of a single peak (ce// = fagg). 
The central region of these clusters do not appear spherical, 
rather this region has some degree of flattening. We see that 
multi-peak systems, mostly bimodal clusters with un-equal 
size sub-clumps (ce// < £ogg), are common for these clusters. 
At low brightness levels, i.e. in the outer regions of clusters, 
the sub-clumps appear in various shapes. In some cases they 
merge forming one system (ce// = eagg), in few cases they 
appear homogeneously distributed (ce// > Cogg), and in few 
cases they form filamentary structure (ce// < Cagg). The 
degree of inhomogeneity (ce// 7 ^ ^agg), generally, is higher 
for X-ray clusters. There is a weak trend that cluster centers 
are more flattened than the outer parts, irrespective of the 
nature of simulation. 


4 RESULTS 

One of the objectives of this paper is to study morphological 
evolution in simulated and optical clusters using M (equa¬ 
tion 4) and e (equations 5 - 7) as quantitative measures. 
These parameters represent the shape characteristics of a 
set of iso-density/intensity contours corresponding to a set of 
density/brightness levels. The levels represent equal interval 
in area, i.e. size in log space, which allow higher resolution 
and hence higher weight to the dense, central region. 

In this study, we emphasize the morphological prop¬ 
erties of the dense, central region of clusters. We analyze 


each cluster at two different threshold levels corresponding 
to radii ^ 0.5 h”^ Mpc and ^ 1 h~^ Mpc where the outer 
radius is within ~ 3 times the core radius (Bahcall 1999). 
For each radius, measurements are relative to the center. 

Cluster images are smoothed by a Gaussian filter with a 
smoothing scale (ss) ^ 50 h~^ kpc. We choose this scale after 
trials with different values. Our experience shows that for a 
scale smaller than ^ 50 h~^ kpc, images contain too much 
noise whereas for a larger scale they become over smoothed. 
We note that smoothing affects the gross morphology with¬ 
out much distortion in the evolutionary trend of the pa¬ 
rameters. The trend holds for both simulated and observed 
cluster samples. 

In Figs. 6 and 7, we show detail properties of adiabatic 
and observed clusters using Me// and tagg. These figures 
show results for ss = 50 h~^ kpc within a radius of 0.5 
h“^ Mpc (panels numbered 1) and 1 h“^ Mpc (panels num¬ 
bered 2). Simulated clusters are shown by (faint) horizontal 
lines at three different redshifts and the optical clusters are 
shown by (dark) crosses. The expressions at each panel rep¬ 
resent the best fit line relating the mean of the parameter 
to the redshift although we note that at each redshift the 
distribution functions are highly non-Gaussian. An interest¬ 
ing feature of these figures is that, at least within 2 < 0.25, 
optical clusters have similar dispersion in both Me// and 
tagg. Similar behavior is also noticed for other parameters 
in this redshift range, irrespective of simulation types. The 
wide spread in multiplicity and projected shape of DM ha¬ 
los and X-ray gas is a clear reflection of different merging 
history (Jing & Suto 2002). Note that the error bar in the 
normalizations (i.e., intercepts) of the best fit lines is less 
than 10% for all parameters (not shown in these figures). 

We quantify the rate of evolution by the slope of the 
best fit line where the rate means either dM/dz or dtjdz. 
For gross morphology, we refer to the normalization of this 
line. We present our final results in Figs. 8 to 13 within ~ 0.5 
h~^ Mpc and ^ 1 h~^ Mpc radii with ss ^ 50 h“^ kpc both 
for simulated and optical samples. For simulations we show 
the best fit line along with its expression in gray color. We 
divide the optical sample into four bins with equal number 
of clusters in each bin. In this case the best fit line is shown 
in dark color. No expression is given for this line. In both 
simulated and observed clusters the error bar represents the 
error in the mean. 

4.1 Comparison among Simulated Cluster 
Samples 

A visual examination of Figs. 8-13 clearly shows an evolu¬ 
tionary trend in cluster morphology since the gross proper¬ 
ties of clusters indeed change with redshifts. We are inter¬ 
ested to determine the significance of this trend of cluster 
properties computed at two different different regions sur¬ 
rounding the cluster center and than to compare it with 
observations. 

Gluster properties within ^ 0.5 h~^ Mpc of the differ¬ 
ent samples are shown in Figs. 8 - 10. In terms of multiplic¬ 
ity, a parameter that probes the number of sub-components 
present in a complex system, we find that cooling samples 
have slightly higher value of multiplicity at all redshifts com¬ 
pared to that in the adiabatic sample. The low abundance 
of single component systems with radiative cooling indicates 
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that the dense, cool core sub-structures are long lived fea¬ 
tures (Motl et al. 2004). We find that the feedback mecha¬ 
nism with cooling makes cosmological systems less clumpy 
than systems without feedback. Energy feedback process, 
most likely, slows down the rate of evolution in the X-ray 
clusters than those in the cooling only samples. However, 
this is quite opposite for the DM halos. Higher multiplic¬ 
ity in the X-ray clusters in the cooling samples indicates a 
possibility of less efficient merging in hot baryonic gas. 

In all simulations multiplicity shows a clear trend with 
redshift: clusters have higher (mean) multiplicity at higher 
redshifts. Cluster multiplicity reflects sub-structures merger 
rate. It decays by the rate at which the cluster can relax, a 
time scale which is roughly equal to the dynamical time. The 
CDM halos host a larger amount of substructure at higher 
redshifts because of lower accretion time as compared to 
the dynamical time (see Zentner et al. 2005). This is the 
reason for the systematic increase in overall sub-structures 
with increasing redshifts. 

In terms of ellipticity, we find that the X-ray clusters, 
in general, are more regular than the halos. The X-ray emit¬ 
ting hot gas is supported by the thermal pressure. Due to its 
isotropic pressure support the X-ray gas becomes homoge¬ 
neously distributed in the background DM potential where 
it evolves. As a result, morphology of the distribution of X- 
ray gas appears more regular. Our results suggest that in 
X-ray clusters the irregular sub-components are distributed 
over a region instead of making a filamentary structure along 
one direction. A comparison of teff with tagg for the halos 
(in all samples) shows that the halo sub-clumps are not dis¬ 
tributed uniformly around the central region. Rather these 
clumps are spread out mostly in one direction forming fila¬ 
mentary structure, as indicated by the larger value of tagg- 

No significant evolution is signaled by Ce// for the DM 
clusters in any of these samples. Recall that this parameter is 
an indicator of shapes of individual components in a cluster. 
Therefore, no evolution means that shapes of isolated com¬ 
ponents in clusters at one redshift appear similar any other 
redshifts. Since it places emphasis on individual component, 
therefore, it is not unusual to find no evolution quantified 
by this parameter. However, shapes of sub-clumps in the 
distributions of X-ray gas change in the cooling simulations 
compared to other simulations. 

Properties of simulated clusters within ~ 1 h“^ Mpc are 
shown in Figs. 11-13. We notice that in this distance, small 
scale structures of simulated clusters do not change signifi¬ 
cantly than what we find in smaller scale. This implies that 
the small sub-clumps can exists up to a Mpc scale and dis¬ 
tributed widely over the cluster body. However, in this scale 
sub-structures evolves a bit faster. We find that individual, 
isolated components become a bit more flattened and their 
evolution is slightly stronger. The overall shape of the clus¬ 
ters, however, is less flattened than the central region and 
evolution is weak in all simulations. 

We take projections along each axis at a time and re¬ 
peat our analysis. Recall that in this case, each sub-sample 
(along each axis) has only 41 clusters. The analysis of these 
sub-samples do not show significant variation from the pri¬ 
mary sample. Therefore, it is unlikely that the overall result 
may have contaminated by the projection effect. We have 
also repeated our analysis using different values for the den¬ 


sity/brightness levels. Apart from a minor change in gross 
morphology, we find similar results for the rate of evolution. 

We summarize our main results as follows: First, the 
DM halos show very similar evolution in all samples of clus¬ 
ters. Second, the X-ray clusters in the adiabatic simula¬ 
tion evolve faster than those with radiative cooling. Third, 
morphology of the central parts of clusters evolve slightly 
strongly than the outer regions. Fourth, feedback processes 
with cooling makes the DM halos slightly more flattened 
and slower in evolution than the cooling only simulations 
(see Kazantzidis et al. 2004 for a similar trend). 

We emphasize that the measured quantities for the DM 
distributions in all three samples are very similar. This is 
a check on the consistency of the simulations and analysis. 
The result is expected as the N-body segment of the sim¬ 
ulations are identical in all three cluster samples with the 
exception of the gas that makes a relatively minor contri¬ 
bution to the total gravitational potential. The LSS of adi¬ 
abatic and cooling clusters are generally similar but their 
small scale structures are determined by the overall clus¬ 
ter properties rather than perturbative interactions (Motl et 
al. 2004). In the adiabatic clusters, the mixing of in-falling 
sub-clumps into the main cluster medium is quicker rela¬ 
tive to the radiative cooling clusters where sub-structures 
can be long lived. This is a reason behind the fast evolu¬ 
tion of adiabatic X-ray clusters. The relaxation time scale 
for collisionless particles is much longer than that of the 
collisional gas particles (Frenk et al. 1999; Valdarnini, Ghiz- 
zardi & Bomometto 1999). Therefore, the DM halos will 
appear not only more elongated than the distributions of X- 
ray gas, but the redshift evolution of their shapes will also 
be slower. More spherical configurations for X-ray clusters 
is also expected from the point of view that intra cluster gas 
is approximately in hydrostatic equilibrium and supported 
by isotropic pressure (Sarazin 1988). The DM, on the other 
hand, appears to be distributed like galaxies as indicated by 
recent observations from gravitational lensing (Fischer et 
al. 1997; Fischer & Tayson 1997; Kochanek 2001; Hoekstra 
2003; Hoekstra et al. 2004) and by high resolution hydro- 
dynamical simulations (Nagai & Kravtsov 2004; Kang et al. 
2005; Maccio et al. 2005). 

There is a consensus based on the observations of X-ray 
clusters that cooling affects the mass distribution apprecia¬ 
bly only in the inner ~ 10% of the virial radius of cluster size 
halos (Sarazin 1986). Contrary to that recent high resolution 
hydrodynamic simulations show something quite interest¬ 
ing (see Kazantidis et al. 2004). Kazantidis et al. show that 
there is a significant difference in overall shape between dis¬ 
sipationless and dissipative simulations which is persistent 
up to the virial radius. The virial mass of their DM halos 
ranges from ~ 10^® to 3 x 10^“^ h”^ Mq which translates 
to the virial radius range ~ 0.26 — 0.82 h“^ Mpc assuming 
Avir{z = 0) ~ 337, h ~ 0.7, and pc ^ 1.87h,^ x 10“^® gm 
cm“® (Kolb & Turner 1990; Zentner et al. 2005). Kazantidis 
et al. present their analysis upto the virial radius. However, 
from the trend seen right at the virial radius, it seems likely 
that it goes a bit further down along the radial direction be¬ 
fore shapes in dissipative and dissipational simulations con¬ 
verge. 

Baryon fraction (fl^ ~ 0.043) in Kazantidis et al. simu¬ 
lations is larger than what has been used in our simulations. 
Therefore, question can be raised whether low baryon den- 
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sity can also produce systematic shift in the shapes of DM 
halos to be robust on scales of Mpc as noted in our work. Our 
simulations use baryon density (^6 = 0.026) and normaliza¬ 
tion of fluctuation spectrum (as = 0.928) which are slightly 
off than the corresponding WMAP values {Qb = 0.044 and 
(T8 = 0.84, Spergel et al. 2003). Baryon density is an impor¬ 
tant cosmological parameter which affects radiative cooling 
and X-ray luminosity at the central region of large virial- 
ized structures. Higher Qb enhances the cooling rate, sub¬ 
sequently making the central region more regular (Sarazin 
1986, Kazantzidis et al. 2004, Springel & White 2004, All¬ 
good et al. 2005, Flores et al. 2005). Recent numerical sim¬ 
ulations show that larger as produces DM halos that are 
more regular in the central regions (Allgood et al. 2005). 
Therefore, we note that cooling is under-emphasized while 
the core DM sub-structure is over-emphasized is our simu¬ 
lations. It may be likely that the offset of Qb and as com¬ 
pared to WMAP would balance each other and our results 
obtained from the cooling simulations would still be repre¬ 
sentative had we been using the WMAP values. 

With cooling only, our simulated clusters of galaxies 
show a large amount of long-lived substructure compared 
to the other simulated samples. While the amount of cool¬ 
ing in this sample is unphysical it represents an interesting, 
theoretical, limiting case. On the scale of the cluster itself, 
the gravitational and dynamical effects of cool, dense cores 
of gas have significantly altered the shape of the clusters to 
length scales comparable to the virial radius (see Figs. 4 & 
5). The perturbation from cool baryonic clumps may thus 
significantly alter model dependent mass maps derived from 
weak lensing studies. The robust substructures present in 
the cooling only sample may also play a role in steepening 
the total cluster mass profile (Maccio et al. 2005) and with 
higher resolution simulations may bound the possible contri¬ 
bution of substructures to strong lensing in clusters. Though 
beyond the scope of the current paper, these connections to 
lensing studies will be pursued in future work. 

Numerical simulations provide interesting information 
on two different aspects of the LSS: 1) shape of the central 
structures in galaxy or cluster size halos, and 2) Change in 
shape of halos with radial distance, irrespective of the nature 
of simulation. Recent high resolution hydrodynamical sim¬ 
ulations has shown quite successfully that hydrodynamical 
phenomena make cluster centers considerably more spherical 
than those in the adiabatic simulations. However, radial de¬ 
pendence of shape is still a controversial issue. While Frenk 
et al. (1988), Bullock (2002), Springel, White & Hernquist 
(2004), Hopkins, Bahcall & Bode (2005), and hydro simu¬ 
lation of Kazantidis et al. (2004) agree that inner part of 
clusters are more spherical than the outer part, the follow¬ 
ing groups of Dubinski & Carlberg (1991), Warren et al. 
(1992), Jung & Suto (2002), Allgood et al. (2005), and hy¬ 
dro simulation of Tissera & Dominguez-Tenreiro (1998) find 
it completely opposite. Our results closely follow the latter 
group. 

Note that radial dependence of shape is not monotonic. 
It changes in a quite complicated way depending on the 
presence of sub-clumps as one can see from Fig. 5. A similar 
trend is also seen in hydro simulations of Kazantzidis et al. 
(2004). 


4.2 Comparison with Optical Clusters 

We have analyzed a sample of AGO clusters within red- 
shift, 2 < 0.25. The sample contains 208 optical clusters 
derived from 10-inch photographic plates taken with the 48- 
inch Palomar Schmidt Telescope (for details of the data ac¬ 
quisition and processing see Trevese et al. 1992; Flin et al. 
1995; Trevese et al. 1997; Flin et al. 2000). Results obtained 
from the optical clusters are shown in Figs. 8-13 using dark 
dashed lines. The summary of our results is as follows: 

• The optical clusters are, in general, more clumpy than 
the simulated DM halos as given by both M^ff and Mmax- 
The sub-structure at the central part of X-ray clusters in the 
RC sample are compatible with the optical clusters, at least, 
within redshift 2 < 0.25. At large radius, the optical clus¬ 
ters include more small scale structures and show stronger 
evolution in sub-structures. 

• The sub-structures of hot baryonic gas evolve much 

strongly in the adiabatic simulation than that in the galaxy 
distribution (dMeff/dz ~ 0.14,0.13 in 0.5 and 1 h“^ Mpc). 
On the other hand, effective sub-clumps of the ha¬ 

los have faster rate in both adiabatic and SFF simulations. 
Feedback process along with radiative cooling make rapid 
evolution in DM halo structures. In terms of Mmax, how¬ 
ever, evolution of the galaxy distribution is always stronger 
compared to all three simulations {dMmax/dz ~ 0.57,0.79 
for 0.5 and 1 h”^ Mpc, respectively). 

• The largest component of the DM halos (probed by 
teff) in the adiabatic and SFF simulations have higher elon¬ 
gation compared to that in the galaxy distribution. In the 
RC simulation we find an opposite trend. The shape of the 
largest sub-clump formed in the distributions of X-ray emit¬ 
ting hot gas in all three simulations are significantly rounder 
than that of the optical clusters (dtaff/dz ~ 0.27, 0.22 in 0.5 
and 1 h“^ Mpc, respectively). 

• The overall shape (probed by tagg) and the strength of 
evolution in optical clusters (dtaggidz ^ 0.28, 0.2 in 0.5 and 
1 h“^ Mpc, respectively) show nice agreement with that of 
the X-ray clusters in dissipative simulations. In dissipation¬ 
less simulation, however, hot gas is systematically less elon¬ 
gated but evolve much strongly than the galaxy distribution 
[diagg/dz ~ 0.28,0.2 in 0.5 and 1 h“^ Mpc, respectively). 

• The shapes of the optical clusters are comparable to 
the halos only in the RC simulation. The halos are slightly 
more flattened and slower in evolutionary process compared 
to the galaxy distributions in the adiabatic case and in sim¬ 
ulation including feedback processes with cooling. 

• The strength of shape evolution given by diagg/dz is 
slightly stronger around the cluster core, in both observed 
and simulated clusters. 

There are several possibilities for optical clusters to be 
more clumpy. First, the choice of smoothing may not be op¬ 
timal for the optical sample. The smoothing scale used in 
our study, therefore, should be taken as the lower limit. Sec¬ 
ond, projection effect due to the background galaxies may 
also play an important role. This effect becomes significant 
as one moves away from the clusters center (Kolokotronis et 
al. 2001). Third, since ACO clusters are selected via richness 
criteria and it has been shown that richness is poorly cor¬ 
related with mass (Girardi et al. 1998; Miller 2004, private 
communication), there is a chance that optical sample may 
be biased toward high mass end of cluster mass range or clus- 
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ters that have gone through recent merger. Massive clusters 
are dynamically less relaxed and hence rich in sub-structure. 
A well defined mass selection criteria needs to be applied for 
a more systematic comparison as recent numerical simula¬ 
tions show that cluster shapes depend on mass, although the 
mass - shape correlation is weak and show large dispersion 
(Bullock 2002; Jing & Suto 2002; Hopkins, Bahcall & Bode 
2005; Allgood et al. 2005). 

Regarding the comparison of the strength of evolution 
we note that the morphological parameters derived for the 
set of simulated halos have uniform statistical weight at all 
redshifts. However, this is certainly not the case for the ob¬ 
served sample as it has considerably more weight toward 
« = 0.0 than the simulated samples (see Fig.s 6 and 7). As 
mentioned earlier the scatter is comparable in both samples 
and, therefore, we believe the choice of binning has less effect 
on the overall outcome of our analysis. 


4.3 Results from Previous Studies on Observed 
Clusters 

In this section we summarize the results obtained from pre¬ 
vious studies on optical and X-ray clusters. Our objective is 
to highlight the fact that different samples of clusters give 
different rates of evolution in cluster morphology. Due to the 
methodological differences, we refrain from making a direct 
comparison with the results of these studies. 

To find cluster shapes, all previous studies follow the 
procedure described in Carter & Metcalf (1980). These stud¬ 
ies, however, differ in adopting weighting factor, threshold 
level, cluster center and smoothing techniques to construct 
galaxy density distribution from spatial distribution. It is 
also important to note that cluster shape quantified by ellip- 
ticity is not uniquely defined. Therefore, to help reader get¬ 
ting a better feelings about the inherent differences of pre¬ 
vious studies we also provide a brief outline of the method¬ 
ologies used in these studies. 

The optical sample of MCM contains 138 ACO clusters 
with 2 < 0.1 which has been compiled from West & Bothun 
(1990); Rhee, van Haarlem & Katgert (1991) and Kolokotro- 
nis et al. (2001). This sample show no significant evolution, 
dtjdz ~ 0.03. 

The former two groups measure cluster shapes from dis¬ 
crete galaxy distribution using method of moments. They 
define the two-dimensional moments as, 

E,,, (®i-a:o)"*(j/.-?/o)" 

Pmn = - - --, (8) 

where xo,yo are the coordinates of the brightest galaxy 
taken as the cluster center, N is total number of galaxies 
within the region which is 3cr above the background noise 
and m, n = 0,1,2. They diagonalize the matrix formed by 
the components /r20,M02, and /in, find the eigenvalues and 
obtain cluster shape using eigenvalues from the relation, 
e = 1 — Ai/Ai, where Ai > A2. Kolokotronis et al. (2001) 
use moment of inertia method for a sample containing 22 
ARM clusters along with their ROSAT counterparts in the 
redshift range, 2 < 0.13. They use Gaussian smoothing on 
galaxy density distribution and define the components of the 


symmetric inertia tensor as, 

-^11 = - xi), I22 = '^WiRi - yi), 

i i 

I12 = /21 = - ^ WiXiyi, (9) 

i 

where Wi is the average cell density within 0.75 h“^ Mpc 
region and R = x^ + y^. After defining inertia tensor, 
Kolokotronis et al. follow similar route to the other groups 
to define shape except that they define ellipticity as, e = 

1 — A 2 /A 1 . 

Plionis (2002) analyze the largest sample of optical 
clusters following the method used in Kolokotronis et al. 
(2001). His sample has 407 APM clusters within a volume 
of 2 < 0.18. The rate of evolution for the Plionis sample is 
dtjdz ~ 0.7. However, if both are combined, replacing the 
common ones by the APM clusters, the rate increases. The 
combined sample of ^ 500 optical clusters with 2 < 0.18 
shows dtjdz ~ 1.06. 

It is rare to find a large sample of X-ray clusters with 
up-to-date ellipticity measurements. The X-ray sample of 
MCM is compiled from Mcmillan, Kowalski & Ulmer (1989; 
hereafter MKU) and Kolokotronis et al. (2001). MKU mea¬ 
sure cluster shape using method of moments from 2D X-ray 
surface brightness images. They adopt the following defini¬ 
tion of the moment, 

y-mn = - - --7-. (10) 

z^ij Jij 

where a;o, yo are the components of the image centroid, xo = 
fij and yo = They determine 

the overall shape of a cluster using the faintest flux level 
available for that object. This sample has 48 clusters with 

2 < 0.1 which is three times smaller than the MCM optical 
sample and an order of magnitude smaller than the APM 
sample. It also has a lower redshift limit than the APM 
sample. The rate of evolution for this sample is dtjdz ~ 1.7. 
The result suggests faster evolution for the X-ray clusters 
than the optical one. Interestingly, a comparison of optical 
and X-ray clusters within Kolokotronis et al. (2001) sample 
show completely opposite trend: galaxy density distributions 
have stronger evolution than the distribution of hot gas. 
The galaxy and X-ray cluster shapes follow a trend where 
flattened gas distribution signals anisotropic distribution of 
galaxies. However, the scatter is large in both relationships. 
It is not clear to us what could be the reasons of possible 
contradictions except the fact that MCM sample is most 
likely contaminated due to different methodologies. Besides 
it is also difficult to make any definite conclusion because 
of the smaller sizes of the samples. A large sample of X-ray 
clusters with better selection criteria and extended to higher 
redshift is needed. 

We reanalyze t — z estimates derived from the APM 
cluster data and the combined sample imposing a redshift 
cutoff 2 < 0.1 in order to be consistent with the redshift 
range of MCM X-ray sample. For these samples, we find the 
rate of evolution as , dtjdz ^ 1.02 and ^ 1.0, respectively. 
We find that evolution of optical clusters accelerates in this 
redshift range but it is still slower than that of the X-ray. 

Flin, Krywult & Biernacka (2004; hereafter FKB) has 
analyzed a sample of 246 ACO clusters for 2 < 0.31. This 
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group use the same definition of moment as in equation 1. 
They use density peak as the cluster center and measure 
shapes at different circular aperture radii ranging from 0.5 
to 1.5 h“^ Mpc with an increment of 0.25 h“^ Mpc. They 
estimate cluster shape at all radii and find no dependence 
of cluster ellipticity on redshift. Interestingly FKB noted a 
decrement of dejdz with radius. They find positive evolu¬ 
tion at radii of 0.5 and 0.75 h”^ Mpc. However, for radii 
> 1 h“^ Mpc, they report negative evolution. The mean of 
their estimates derived from these five radii shows e « 0.22 
and di/dz ^ 0.013. For 2 < 0.1, their result also indicates 
weak evolution. We use this sample of optical clusters for 
our analysis (see §4.2) but with a reduced number (208) of 
clusters. The reduction is made after visual inspection and 
it is due to the removal of clusters images that appear either 
small or close to the boundary. 

It should be noted that MCM and APM samples em¬ 
phasize cluster morphology in two different regions. The 
MCM sample excludes any study with radius less than 1 
h”*^ Mpc and includes the estimate of ellipticity within ~1 - 
2 h“^ Mpc from the cluster center. The APM sample, how¬ 
ever, provide information on cluster shape within 0.75 h~^ 
Mpc of the center. Therefore, care must be exercised in in¬ 
terpreting and comparing results of observed clusters with 
simulations if both are not analyzed under the same mea¬ 
surement technique. Unfortunately the studies of Floor et 
al. (2003) and FMM (2004) has ignored this fact. 

In spite of differences in the evolution of cluster mor¬ 
phology, optical samples are consistent with one another 
atleast in one case: shape of galaxy density distributions 
evolve strongly in the central region (Plionis, 2002) than 
that in the outer part (MCM, 2001). Interestingly our re¬ 
sults are also consistent with this trend. For X-ray clusters 
this trend has yet to establish. 


5 CONCLUSIONS 

Numerical simulations provide an unique opportunity to fol¬ 
low the hierarchical nature of the LSS formation in both lin¬ 
ear and nonlinear regimes (Frenk et al. 1985, 1988; Quinn, 
Salmon & Zurek 1986; Efstathiou et al. 1988). In order to be 
representative of the reality, results from simulations should 
agree with observations. Observations provide evidence of 
morphological evolution in galaxy-clusters (Melott, Cham¬ 
bers & Miller 2001; Plionis 2002; Jeltema et al. 2005), sim¬ 
ulations should show similar trend. Besides, in the CDM 
model luminous galaxies are associated with the DM sub¬ 
halos which reside in bigger parent halos, closely associ¬ 
ated with galaxy clusters. According to this model statistical 
properties of galaxies, e.g. mass, sub-structure, shape etc., 
would show a similar trend to that of the sub-halos while 
X-ray emitting hot gas would have different properties than 
galaxies and sub-halos. A statistical analysis of various prop¬ 
erties of halos, galaxy clusters, and X-ray gas could provide 
clues to find possible biasing of luminous galaxies toward 
DM sub-halos and whether or not they have any correspon¬ 
dence with the distribution of hot gas. With this in mind, we 
have studied redshift evolution of cluster morphology sim¬ 
ulated, respectively, in the adiabatic limit, with radiative 
cooling, and with star formation including SN feedback at 
three different redshifts, 2 = 0.0,0.10, and 0.25. For com- 


parision we have also studied a sample of observed clusters 
containing 208 ACO clusters within redshift, 2 < 0.25. 

Since observed clusters are projected along the line of 
sight and lack the full three dimensional information we, 
therefore, use projected simulated clusters. Each cluster im¬ 
age is a 8 h~^ Mpc frame containing 360 x 360 pixels. Clus¬ 
ters are analyzed at two different density/brightness thresh¬ 
old levels corresponding to radii 0.5 and 1 h~^ Mpc from 
the cluster center. To quantify morphological evolution we 
use multiplicity and ellipticity as two different probes that 
are sensitive to cluster sub-structures and shape. 

Our results indicate that optical clusters have, in gen¬ 
eral, more sub-structures than simulated halos and X-ray 
brightness distributions. Cluster components, in both ob¬ 
served and simulated clusters, evolve with redshifts and the 
evolution is different at different regions from cluster centers. 
In terms of total multiplicity (Mmax), observed clusters have 
stronger evolution compared to DM halos. The X-ray bright¬ 
ness distributions, however, show steeper evolution (than 
that of galaxy clusters) in dissipationless simulation. 

We find that in terms of overall shape, simulations do 
model the observed universe in an interesting way. The sim¬ 
ulated clusters evolve with redshift, consistent with the hier¬ 
archical formation scenario. However, observed clusters ap¬ 
pear to be slightly more flattened at higher redshift than 
the simulated one indicating slower evolution in simulated 
objects. This may reflect some form of incompleteness in 
our understanding in simulating the LSS. Our results differ 
from those of FMM (2004) who reported that the evolu¬ 
tion in the simulated cluster shape is significantly slower 
than the observed one. We not only find stronger structural 
evolution in simulated clusters, but also find that observed 
cluster shapes appear to be consistent with dissipative sim¬ 
ulations, at least, in the redshift range 2 < 0.1. The discrep¬ 
ancies noted in FMM is due to the different redshift ranged 
probed as well as intrinsic methodological differences while 
comparing simulations with observations. 

We note that on one hand shapes of optical clusters 
seems to be compatible with both the halos and X-ray 
brightness distributions, one the other hand, both of these 
components appear to be less clumpy than the distribution 
of galaxies. Therefore, it seems puzzling whether or not there 
is any correspondence between the DM halos and galaxies. 
The existence of any such correspondence is still a matter 
of ongoing debate as there are conflicting results based on 
systematics of numerical simulations such as nature of simu¬ 
lations (dissipationless or dissipative) and the effect of mass 
and force resolution (see Maccio et al. 2005). In the context 
of the CDM model we would expect that the optical clus¬ 
ters would have similar morphology and evolutionary trend 
to that of the halos and would be different than the proper¬ 
ties of the distribution of hot gas traced in the X-ray region 
of the spectrum. 

Within the uncertainties and systematics involved in 
our optical sample, the results indicate that the properties 
of optical clusters do not exactly represent either the distri¬ 
bution of the halos or that of the X-ray emitting gas in any of 
the simulations. We find offsets in the measured parameters, 
such as multiplicities and ellipticities, between observations 
and simulations, and are unable to find any clear signature 
of DM-galaxy biasing based on our morphological analysis. 
This may be an indication, although in no way conclusive. 
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of the fact that these components of the LSS may repre¬ 
sent intrinsically different populations, and galaxies may not 
trace the DM distributions (see Gao et al. 2004a,b; Nagai & 
Kravtsov 2005). However, this is merely a speculation and 
we stress that care must be exercised in interpreting our 
results as one must be careful in selecting proper measures, 
radius, mass range, and most importantly, well defined sam¬ 
ples of clusters to have unbiased and meaningful results in 
any morphological analysis comparing observations and sim¬ 
ulations. 

We hnd that the measurements from different samples 
do not agree on the evolution rate. Take, for example, op¬ 
tical clusters with 2; < 0.1, and radius, 0.75 h~^ Mpc. In 
this case, the APM sample shows dejdz ~ 1.02. FKB, on 
the other hand, finds much weaker evolution, dtjdz ^ 0.2. 
As mentioned in FKM, the discrepancy may be due to dif¬ 
ferences in adopting cluster centers, smoothing, and applied 
method of shape determination. 

A preliminary analysis of a sample of 800 clusters con¬ 
structed from the Sloan Digital Sky Survey (SDSS) shows 
that ellipticity evolution of optical clusters, for 2 < 0.1 
and within ~ 1 h~^ Mpc, is weaker than that of the 
APM clusters. The result indicates that clusters with dif¬ 
ferent mass limits evolves differently. Large, massive clusters 
(M ~ 10^® M0) have stronger evolution compared to the less 
massive clusters (M ~ 10^® — 10^^Mq) (C. Miller 2004, pri¬ 
vate communication). This is an interesting observation. If 
it is confirmed then the scaling relation between axes ratios 
and mass noted in simulations (Bullock 2002; Jing & Suto 
2002 ) must be modihed to be consistent with observations. 

The SDSS sample is uniform with a well documented 
selection function and high degree of completeness. We may 
then infer that the cluster samples discussed previously 
have less uniformity in mass range: the APM catalog and 
FKB samples are biased toward massive clusters whereas 
the MCM samples contain more less massive clusters. The 
discrepancy may also arise from the techniques applied in el¬ 
lipticity estimates (see also Flores et al. 2005 in this regard). 
Unfortunately we are unable to check the evolution strength 
- mass relation for our optical sample because, apart from 
an approximate range, no well defined criteria has been used 
to sort clusters into different mass bins. 

The discrepancy in the optical samples is an indication 
of different selection criteria used to construct the catalogs. 
Larger and more complete catalogs obtained from the SDSS 
and XMM-Newton survey may be able to shed more light 
into this issue. It is also likely that numerical simulations 
may lack crucial physics that needs to be included (see FMM 
for discussion). In the future we will analyze clusters simu¬ 
lated with various gas physics, e.g. thermal conduction and 
AGN heating, and compare them with the SDSS clusters. 
The results of these studies may give us some clues to gain 
better insight of the current discrepancy. 
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Figure 1. Contour plots of toy clusters at different brightness 
levels (in arbitrary scales). The multi-modal clusters have clumps 
with different peak brightness. For all clusters the outer line rep¬ 
resents the percolation level where the sub-structures merge and 
form a single, large system. 
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in the text. The x-coordinates of these legends are chosen only for 
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highest peak along the x-axis corresponds to the Mmax whereas 
Meff is obtained after averaging along the x-axis. See text for 
details. 
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Figure 4. Multiplicity (M) as a function of contour area (A 5 ) for 
a selection of clusters at 2 = 0.50,0.25, 0.10, and 0.0. Two clus¬ 
ters from each redshift are shown. Dark, gray, and faint solid lines 
represent respectively, the adiabatic, radiative cooling (RC), and 
star formation with feedback (SFF) samples. The dark matter 
(DM) and X-ray clusters are shown on the left and right panels, 
respectively. Multiplicity is, in general, greater than 1 in the en¬ 
tire redshift range for clusters simulated with RC (medium line) 
indicating a slower evolution than in the adiabatic sample (dark 
line). Redshift 2 = 0.5 is taken for demonstration purpose only. 
See text for details. 
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4. Solid and dotted lines are used to represent eagg and Cg//, 
respectively. The color style is similar to Fig. 4. In most cases 
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peak (i.e. Cg^y = ^.agg) whereas in the outer regions sub-clumps 
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trend is weak. See text for details. 
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Figure 6. A detailed comparison of the estimate of Meff for 
the adaibatic DM (left panels) and X-ray (right panels) clusters 
and the optical sample with ss=50 h~^ kpc within 0.5 h~^ Mpc 
(panels 1) and 1.0 h~^ Mpc (panels 2) radius. Simulated clusters 
are shown by (faint) horizontal lines at 2 = 0.25,0.10,0.0 and 
the optical clusters are shown by (dark) crosses. The expressions 
represent the best fit lines for observation (dark; top line in each 
panel) and simualtion (faint; second line in each panel). 
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Figure 7. A detailed comparison of the estimate of eagg obtained 
from the adaibatic DM (left panels) and X-ray (right panels) clus¬ 
ters and from the optical sample. Presentation style is similar to 
Fig. 6. 
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Figure 8. Adiabatic sample with 50 h~^ kpc smoothing (ss) 
within 0.5 h”^ Mpc radius. Dark and gray lines are used for opti¬ 
cal and simulated clusters, respectively. The error bar represents 
the error in the mean. The expression at each panel relate the 
evolution of the mean value of the parameter of with redshift. 
The strength of evolution for optical clusters are: dM^ff ~ 
0.14, dMmaxIdz 0.57, de^jjjdz ~ 0.27, and deaggjdz 0.28. 
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Figure 9. Radiative cooling (RC) sample with 50 kpc Figure 11. Adiabatic sample with 50 h~^ kpc smoothing (ss) 

smoothing (ss) within 0.5 h~^ Mpc radius. Presentation style is within 1.0 h~^ Mpc radius. Presentation style is similar to Fig. 

similar to Fig. 8. 8. The strength of evolution for optical clusters are: dM^jj!dz ^ 

0.13, dMmaxIdz 0.79, deeffjdz ~ 0.22, and deaggjdz 0.20. 
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Figure 10. Star formation with feedback (SFF) sample with 50 

kpc smoothing (ss) at 0.5 h~l Mpc radius. Presentation style Figure 12. Radiative cooling (RC) sample with 50 h~l kpc 
is similar to Fig. 8. smoothing (ss) within 1.0 h”^ Mpc radius. Presentation style is 

similar to Fig. 8. 
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Figure 13. Star formation with feedback (SFF) sample with 50 
h”^ kpc smoothing (ss) within 1.0 h~^ Mpc radius. Presentation 
style is similar to Fig. 8. 
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